arXiv: 1508.01445v2 [hep-ph] 22 Sep 2015 


Factorization for radiative heavy quarkonium decays into scalar 

Glueball 


Ruilin Zhu * 1,2,3 

1 INPAC, Shanghai Key Laboratory for Particle Physics and Cosmology, 
Department of Physics and Astronomy, 

Shanghai Jiao Tong University, Shanghai 200240, China 
2 State Key Laboratory of Theoretical Physics, Institute of Theoretical Physics, 
Chinese Academy of Sciences, Beijing 100190, China 
3 CAS Center for Excellence in Particle Physics, Institute of High Energy Physics, 
Chinese Academy of Sciences, Beijing 100049, China 

Abstract 

We establish the factorization formula for scalar Glueball production through radiative decays 
of vector states of heavy quarkonia, e.g. J/tf, if(2S) and T (nS), where the Glueball mass is 
much less than the parent heavy quarkonium mass. The factorization is demonstrated explicitly 
at one-loop level through the next-to-leading order (NLO) corrections to the hard kernel, the non- 
relativistic QCD (NRQCD) long-distance matrix elements (LDMEs) of the heavy quarkonium, 
and the light-cone distribution amplitude (LCDA) of scalar Glueball. The factorization provides 
a comprehensive theoretical approach to investigate Glueball production in the radiative decays 
of vector states of heavy quarkonia and determine the physic nature of Glueball. We discuss the 
scale evolution equation of LCDA for scalar Glueball. In the end, we extract the value of the decay 
constant of Scalar Glueball from Lattice QCD calculation and analyze the mixing effect among 
/o(1370), / o (1500) and / o (1710). 
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I. INTRODUCTION 


Glueball, as a kind of color-confined state with two or more than two gluons, is one of 
the most important expectation in Quantum Chromodynamics (QCD). The non-Abelian 
interactions among gluons tell us the existence of Glueball and there is no hidden principle 
to forbid this kind of state up to now. Many theoretical approaches have been employed 
to investigate one of the most important quantum numbers of Glueball, i.e. its mass, e.g. 
Lattice QCD pi, QCD Sum Rules [5HE], Supergravity Dual [9], Top-down Holographic 
Dual [ffij, Rotating Closed Strings [H|, MIT Bag Model [12]. In these approaches, the 
scalar Glueball is expected generally to populate the low enery region from lGeV to 2.2GeV, 
which is also a region rich in qq states. A famous criterion to distinguish Glueball from the 
conventional qq states is that the width of Glueball is narrow from the large- N c argument [EH 
ns], where the decay width of the qq states is proportional to 1 /N c while the width of Glueball 
is proportional to 1/1\T?. Another criterion is that Glueball with non-zero spin (J 7 ^ 0) is 
blind to quark flavor, while scalar Glueball with the quantum number 0 ++ decays to qq is 
chiral suppressed pa nz] and thus its decay to ss is favored than to uu or dd. 

The search of Glueball can be carried out in many experiments: pp collision, nN scattering 
on polarized/unpolarized targets, central hadronic production, D-meson decay, and radiative 
decays of vector states of heavy quarkonia, i.e. J/tp , ip(2S) and T (nS) [T5HZT] . A great deal 
of data has been obtained and it is believed that a large possibility of Glueball component in 
the scalar state / 0 (1500) or / 0 (1710) [IB]. But a pure Glueball where only gluons contribute 
the constituent has not been observed or verified up to now. The study of glue-rich processes 
shall be greatly helpful to hunt the signal of Glueball. The radiative decay of vector states 
of heavy quarkonia through V{1 ) —y G + 7 is one of the most important processes which 

shall provide a platform to systematically investigate the properties of Glueball. 

In the case of the mass squared of produced Glueball greatly less than that of the parent 
heavy quarkonium, i.e. mP G « rriy, a large momentum is transferred, the final Glueball and 
photon both run in the light-cone, and the light-cone factorization can be well-employed. 
The soft and collinear physics is only contained in LCDA of Glueball and the NRQCD long¬ 
distance matrix elements of the parent heavy quarkonium. The light-cone operators populate 
the collinear subgroup of four-dimensional conformal symmetry, and the LCDAs of Glueball 
which incorporate the gluon distribution with certain light-cone momentum fraction can 
be defined accordingly. The scale dependence of the corresponding physical observable is 
governed by evolution equation. The evolution equation of LCDAs can be understood as 
the renormalization group equations for the light-cone operators. Some pioneer works on 
the evolution equation of distribution amplitudes for exclusive reactions at large momentum 
transfer can be found in Refs. |22l [27] . The study of Glueball from radiative Upsilon decay 
based on soft-collincar effective theory (SCET) can be found in Ref. [25]. 

In the paper, we establish the light-cone factorization for scalar Glueball production in 
the radiative J/ip decay, which will be tested at next-to-leading order. The factorization 
formula is also valid in other vector heavy quarkonium decays, e.g. ip(nS) and T (nS). 
Through the establishment of the factorization formulae, it is conveniently to perform a 
systematical phenomenology analysis and open a new and clear widow to investigate the 
properties of Glueball. 

We first give the definition of leading-order LCDAs for scalar Glueball. Since the Glueball 
and the flavor-singlet qq state have identical conformal spin, they will mix each other by 
renormalization, which is ananogous to the quark-gluon splitting behavior in Proton. Thus 
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we define a two dimensional light-cone distribution amplitude 


$(u) = 



( 1 ) 


where the twist -2 LCDAs of <j)i(u ) can be written as follows in terms of quark and gluon 
fields 


(p q (u) 


g*( 2u_ l) fc+z / 2 

27T N q 


(G(k)\%(-z-/2)L ij (-z-/2,z-/2)^ j (z-/2)\0), 


( 2 ) 


<t>g{u) = 


g i(2u—l)k + z /2 

27t N g u{ 1 — u) 


g^{G(k)\G^(-z-/2)L ab (-z~ /2,z~ /2)\0),(Z) 


where the resummation of all order soft and collinear gluon radiation from quark or gluon 
field is summarized into the related gauge link, which also ensure the gauge invariant of 
the defined matrix elements. For a vector p, the light-cone component is given by p M = 
(p + , p~ , p 1 , p 2 ) with p + = (p° +p 3 )/-s/2 and p~ = (p° — p 3 )/V 2. The u defined above is 
the momentum fraction in plus direction for one gluon in Glueball. The factor iVj satisfies 
N q — 1 and N g = k + . We define two light-cone vectors n and n with n M = (1,0, 0,0) and 
n^ = (0,1, 0,0) in the light-cone frame. The tensor factor (f^ can be written as (/“' = 
For the gauge link, we have 


L( Xj y\ — p e 'sJ 8 1|i! M(‘ # (H s +!f) ) 


( 4 ) 


where A^ = A^T a in the fundamental representation for the links between quark and anti¬ 
quark; while A M = if abc A^ in the adjoint representation for the links between gluon and 
gluon. 


II. FACTORIZATION FORMULAE 


We consider the radiative decay of J/ip in its rest-frame [29TF3T] 


J/V>(P)-►<?(*:)+ 7 (g), 


( 5 ) 


where the related momenta are given in the brackets. The momentum of Glueball can be 
written explicitly as (0(k + ), 0(k~), O(Aqcd), O(Aqcd))- Since the charm quark is heavy, 
the produced Glueball has a momentum of order of m c , and satisfies 


kr 

k+ 


m 


G 


m 


JN> 


0.30, 


Aqcd V^Aqcd n 1 . 

——— = - ~ IJ. 14 

k + mj 


( 6 ) 


where we assume that me ~ 1.7GeV, A q C d ~ 0.3GeV. If the parent heavy quarkonium 
becomes to T, the corresponding ratios of k~/k + and Aqcd/^ + dramatically decrease to 
0.03 and 0.04 respectively. Thus the light-cone factorization is more practical in the T 
radiative decay. 
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We only consider the QCD corrections here, so the S-matrix element for the decay is 


(t(?)G(A-)|S| J/V-) = -ieQ c e'“(q) J d i xe il ‘ x {G(k)\c(x) 1 ,c{x)\.J^(P )), (7) 


where Q c is the electric charge of charm quark, c(x) is the Dirac held for the charm quark, 
is the polarization vector for the photon. At the leading-order twist defined in Eq. (J3]) , at 
least two gluons bound to Glueball. The corresponding contribution to the S-matrix element 
is 

^(q)G(K)\S\J/^) = -i^ e Q c g 2 s e* a {q) J d*xd*yd*ze* x 

x {G{k)\T[c{x) r ) tl c{x)c{y) r ) ■ G(y)c{y)c{z )7 • G{z)c{z))\J/^{P)) ,( 8 ) 


We can calculate the T-ordered operator product by Wick-contraction and we use the ex¬ 
pansion of the heavy quark relative velocity v [32j 

(OlaMcjMI J/V) = ~l {P + litP~)ji ( 0 | X V^| + 0 ( 1 ?), ( 9 ) 

where i/j and x^ are the NRQCD operators to annihilate the quark and anti-quark respec¬ 
tively. Since the heavy quark relative velocity squared is around 0.3 for J/-0 and 0.1 for T|33j. 
we neglect the contribution from higher orders of v 2 in this paper. By the simplification, 
the amplitude can be written as 

(^(q)G{K)\S\J/i)) = i^eQcg 2 £* a (q)(27r) 4 5(P - k - q)(0\x ] cr l3 'il?\J/'ild} 

x J ^r^(k,q 1 )M aPtl „(P,k,q 1 ), (10) 

where 

^ v (k,Qi)= I d 4 xe-^ 1 '^^- 91 > y (G(A;)|G' a ^(^)G a ’ y (y)| 0 ), ( 11 ) 

the function T MI/ (fc,gi) incorporates the non-local interactions among two gluons and Glue- 
ball, while the function M a g IJM (P, k, q\) is a perturbative kernel, which can be calculated 
order by order. At twist -2 level the function T AtI '(fc,g 1 ) can be simplified into 

= (2n)‘5(q;)S 2 ( qil _) . 1 -r afFoM . ( 12 ) 

u{u — 1 ) 


Fo(u) = 


2n (k 


+D 


dx-e- i(1 - 2u)t *‘ (G(k)\Gl (-x-)G%(x-)\0). 


(13) 


The naive factorization in Eq. (10) is valid in tree-level. At NLO and beyond NLO, 
the factorization should be corrected to including the fluctuation between the gluonium 
component and the qq flavor-singlet component. 

Leaving all possible Lorentz invariant construction, the amplitude of J/-0 radiative decays 
to scalar Glueball can be factorized into 


iM 


ieQ c g 2 s (0\x ] (T^\J/^) 

x du dt(ej/^ ■ £^m 2 c H 0 (u,v, y) + £j/^ ■ qe^ ■ PH x (u,v, n))&(t, /i), 

Jo Jo 

( 14 ) 
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with 


(OlxWI-W) = ^J/i,(v,n){0\x , rrip\J/ip) r , 


(15) 

(16) 


Note that (0\x^cr'^\J/'^) r here is the matrix element after renormalization, which is isolated 
to the renormalization of LCDAs for Glueball. Tj/^ is the renormalization factor of LDME 
for J/i]}. r is the renormalization factor with 2x2 matrix elements, which can be calculated 
through the renormalization of LCDAs of Glueball. Ej/^ and e 7 are the polarization vectors 
of J/-0 and the radiated photon, respectively. Hi is the hard kernel with two components. 
The factorization formula can also employed to T —» y + G by replacing Q c —> Qb, m c —> mb, 
£j/y -t e r , and m J/v , ->• m T . 

The typical Feynman diagrams at both tree and one-loop level contributing to the hard 
kernel for a heavy quarkonium radiative decays to Glueball can be found in Fig. [TJ Other 73 
diagrams for gg'y final states can be obtained by exchanging the outgoing gluons or inverting 
the quark line in Fig. [lj For tree level, other 3 symmetrical diagrams can be obtained by 
inverting the direction of the quark in the first line. For one-loop level, there are another 5, 
3, and 1 pentagon diagrams respectively compared with the typical diagrams in the second 
line. There are another 11, 7, 3, and 0 box diagrams respectively in the third line. And 
there are another 17, 11, and 1 triangle diagrams, another 11 self-energy diagrams in the 
fourth line. While the diagrams in the fifth line denotes the contribution to qq'y final states, 
and another 3 box diagrams and 1 pentagon diagrams are not shown. 
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FIG. 1: Typical Feynman diagrams for the hard kernel of a heavy quarkonium radiative decays to 
Glueball. 
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After the calculation, we can easily get the LO hard kernel 


flf = I 0, 




N c ml^fnJJ^u(l - u ) 


H S 0) = I 0 , 


N c ml v / 2mjJ^uu 


( 17 ) 


The NLO hard kernel can also be obtained after considering the held renormalization and 
counter-term. It can be written as 


H 


(i) 


y/2, 


CXc 


f rif(2u — l)(nlnu + ulnu) 


e 2nN c mlyfmjJ^uu 


{ 


2 uu 


-en f HS, 


Ca(2u 2 — 2u + l)(ulnn + ulnu) /3 0 C F . 2 

- - 1 ^ 1 A — ^ — 1 7T ( —- 

uu 2 Av 2m c v 


I 1 \2e 




0 1 


H W = —-H 


(1) l „ „ 

0 j 


(18) 


where the coefficients Hq, Hq which contribute to the finite term is presented in the Ap¬ 
pendix. 

The renormalization factor T ju can be obtained by the renormalization of the naive 
LDMEs (O|xl<7"0| J/ift). At one-loop level, it can be written as [33] 


T®, = 1 

J/i> 


(i) a s C F ( 2 . ,1 , ft 

1 \!.l = —- I vr — m {—b In 


4vru 


Amlv 2 


(19) 


The factor T can also be obtained by the renormalization of LCDAs for Glueball. Cal¬ 
culating the LCDAs defined in Eq. (JT|), we can easily get the tree-level result 


r ( 0 ) (M,h) = 


5(u — t) 0 
0 5(u — t ) 


( 20 ) 


At one-loop level, the related Feynman diagrams can be found in Fig. [2j and we have 

( 21 ) 


(i) c (i) 

99 


r(U( u t n) = ^(ii) 2e i ( Sqq s> 

1 27t hq e 1 - 


cd) cC 1 ) 

J 93 d 99 


with 


S^(u,t) = Cp J ^1 + Q(t - u) + (u -b U,t -A t) , 


Sql\u,t) = 2n f T F ^{2u-t- l)0(t-u) - {u —y u, t —y t ) 

L L 


S^(u, t ) = C F — (2 1 — u) 9{t — u) — (u —y u,t —y t) , 
S£ } (M) = ^ ^ + 2 (u + t(1 + 2u))^ 0(t - u) 

+^6(u — t) + (u -y u, f -> t) , 


( 22 ) 
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FIG. 2: The renormalization of LCDAs for Glueball at one-loop level. 


where {3 0 = HN c /3 — 2nf/3, the group factors N c = Ca = 3, C F = 4/3 and T F = 1/2 for 
SU C (3), and the plus function is defined as 

F(x,y) + = F(x, y) — S(x — y) [ dzF(z,y). (23) 

Jo 

Note that the evolution kernel of LCDA is silimlar to that of the corresponding non-local 
operators, the latter ones have been investigated in Refs. [261 EH Eg. 


III. EVOLUTION EQUATION FOR LCDAS OF SCALAR GLUEBALL 


The study of the universality of LCDAs for Glueball which describes the long-distance 
interactions effects is a crucial and also interesting issue [361 E7]. One can find that all 
the divergences are cancelled out between the hard kernel and the renormalization factors 
for both LCDAs and NRQCD LDMEs. We now turn to the scale evolution of LCDAs for 
Glueball. The reason is that logarithms of the form (a s \i\(m 2 cb /ijfy) n , where /i 0 ~ 1 GeV 
denotes the scale at which nonperturbative physics of the LCDAs exists, are large and must 
be resummed to all orders. By the renormalization equation, we can resum these large 
logarithms. The scale dependent equation of LCDA for Glueball reads 

/J 2 ) = V(u, t, a s {fi 2 )) <g) $(t, /i 2 ), (24) 

where the evolution kernel V is 


V = -T -1 ® 



Ois (P 2 ) 
2n 


cA) 

Oqq 

cl 1 ) 


cl 1 ) 

cl 1 ) 

°99 


+ ^(«D • 


(25) 
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The eigenfunctions of Eq. (24) after the renormalization are the Gegenbauer polynomials 
and hence the LCDAs of Glueball <fr(w,/i 2 ) possess the expansion J28J j38] 


®q(u,V 2 ) = 6u(l-u)fg <(^ 2 ) C 'n /2 ( 2M - !)) 


n=l,3,... 


*»('“. = 30« 2 (1 -uf! a ( 1+ ^2 <(m 2 )CS-i(2u - 1) ) , 


(26) 


n=3,5,... 


where we omit the even n series as a consequence of the symmetry of LCDA for scalar 
Glueball, i.e. $ 3 (-u, /i 2 ) = $ s (l — u, /i 2 ) and /j, 2 ) = — 40(1 — w,/i 2 ). The Gegenbauer 
momentum a n also obey the renormalization group equation 




0 a s (fi 2 ) 7 V 1? 


ry99 ry99 
In in 


«A 2 ) 


d/g 2 \ a 9 n (n 2 ) ) 27r 

where the anomalous dimensions that govern the evolution of the LCDAs are 

- 4 tj)(n + 2) - 4V>(1) ) , 


(27) 


y)? 9 — Cp I 3 + -— 

,n V (n + l)(n + 2) 

qg _ 24n/T F (?7 2 + 3n + 4) 

n(n + 1) (n + 2) (n + 3) ’ 
C F (n 2 + 3n + 4) 

3(n + l)(ra + 2) ’ 


r^gq — 
In 


In = C A ( -4 i>(n + 2) + 4-0(1) + 


0o 8 (n 2 + 3n + 3) 

2C0 n(n + l)(n + 2)(n + 3) 


(28) 


where 0(a;) is the digamma function. 

After solving the renormalizaiton group equation, the resummed Gegenbauer momenta 
read 


«A 2 ) = 




gq 

n 


<(j 4) {K-I 9 n 9 ) 


OtJjr) 

.<Xs(l4 


2)1 2A+//3 0 


<(/i 2 ) = g (<40o) 


a s (/r 2 ) 


where are the eigenvalues 


2A+//3 0 


°n (Mo) 


«s(/io) 


- «n (>o) (A n - 7n®) 
2\ 1 2A^//3 o' 


_«s(/io) 


2\ 1 2A„ //3q n 


A± = 2 ( 7f + 7^±^) , 


with 5 = \/(yn 9 — yn 9 ) 2 + 47n 9 7n 9 - And a± are the eigenvectors, with 

(A 2 ) = <(^ 2 ) 7 n 9 - <(/^ 2 )(An - 7 f) • 


(29) 


(30) 


( 31 ) 


We show the LCDAs of scalar Glueball in Fig. [3j where we take af(/r 0 ) = 0.2 and af(/xo) 
0 with /i Q = lGeV as input, and evolute it into another scale. 

























FIG. 3: The light-cone distribution amplitude for scalar Glueball, where we only consider the first 
Gegenbauer momentum and resum the corresponding large logarithms in a|(/i) and a\(n), using 
H o = lGeV and af(/ro) = 0.2, a\{n o) = 0 as input. The asymptotic form of <j> g (u) is renormalized 
to 30 u 2 (1 — u) 2 . 

IV. PHENOMENOLOGICAL DISCUSSIONS 

In this section, we will employ the above factorization formulae and analyze the phe¬ 
nomenological results confronting the recent BESIII and CLEO data [20| [2T]. The branching 
ratio of J/ip —» G T 7 can be written as 

m 2 T/ , — m 2 r 

BU/^^G + j) = -> G + t)| 2 . (32) 

16vrr j/^m 3 ^ 

At first, we note that the decay width of J/-0 —» G + 7 has been studied by the CLQCD 
Collaboration within the framework of quenched Lattice QCD |4T] , which gives 

-A G + 7 ) = (3.8 ± 0.9) x 10 " 3 . (33) 

In the following we will adopt the value of parameters from PDG2014: [42J rrtj = 
3.0969GeV, T j/^ = 92.9keV, m^s) — 3.686GeV, E^s) = 286keV, rriy = 9.4603GeV, 
Tx = 54.02keV, my^s) = 10.023GeV, Txjgs) = 31.98keV, mx( 3 S) = 10.355GeV, Tx( 35 ) = 
20.32keV. The values of LDMEs for heavy quarkonia are extracted from their electric widths 
at NLO as Ref. [43], which reads as (O|xtcr0| J/ip) = 0.6408(GeU) 3 / 2 , (O|y t cr'0|'0(2S')) = 
0.4975(GeU) 3 / 2 , (0|x f cr^|T) = 1.710(GeU) 3 / 2 , (0|xW|T(2S)) = 1.2502(GeU) 3 / 2 and 
(O|x' l 'cr'0|T(35')) = 1.099 (GeV) 3 / 2 . The heavy quark mass is adopted as m c = 1.5GeV and 
nib — 4.8GeV |44[|45jJ. From Fig. [3j one can see that the 3> g («) is small, so we will ignore its 
contribution in the following. We take the mass of scalar Glueball as mg = 1.710GeV from 
Lattice QCD [3], then we can extract the decay constant of scalar Glueball, which reads as 

f g = 0.0386l°;^ 9 7 GeV, (34) 
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where the uncertainty is from both the Lattice QCD and the running coupling constant. 
Note that the result is a little smaller than the prediction (0.10 — 0.13)GeV from QCD sum 
rule [3Uj. We will get a more precise and reliable result for the decay constant of Glueball, 
if we have the Lattice result for Y(nS') —» G + 7 , where the NRQCD(^)LCDA factorization 
becomes more solid. We can also predict the branching ratios of scalar Glueball from other 
vector heavy quarkonia. They are 


B(ip(2S) G + 7 ) = (5.9l 3 ; 4 ) x 10 “ 4 , 

£(T^G + 7 ) = (1.3l°;D x 10" 4 , 

B(T(2S)^G + 1 ) = (1.01°;°) x 10“ 4 , 

B(T(3S) G + 7 ) = (1.2l°;£) x 10 “ 4 . (35) 


Next we will consider the mixing among scalar Glueball and scalar qq states. There are 
many scalar mesons with masses lower than 2GeV, which can be classified into two nonets: 
one nonet with mass below lGeV includes / 0 ( 500), / 0 ( 800), Jig (800) and a 0 (980); the other 
nonet with mass above lGeV includes iFg(1430), a 0 (1450) and two scalar mesons [131 SSI- 
One can see that not all three isosinglet scalars /o(1370), /o(1500) and /o(1710) can be 
accommodated in the qq nonet picture at the same time. One of them can have a large 
possibility of Glueball component. Denoting nn = (uu + dd)j\J 2, we write the mixing 
formula [39l @0] 

|/o) = OLi\nn) + Pi\ss) + Pi\G), (36) 

where /g can be one of / 0 (1370), /o(1500) and /o(1710). 

According to the lattice calculations m, Lee and Wcingarten found that /o(1710) is 
composed mainly of scalar Glueball. It is reasonable when we take the current exper¬ 
imental data into account. From PDG2014, we find that B{i]){J/il)) —» / 0 (1710)7) is 
a large value, which will be 1.56 x 10 -3 when considering B(ip(J/ip) —> / 0 (1710)7 — > 
7 T 7 T 7 ) = (4.0 ± 1.0) x 10 -4 , B(j/j(J/ip) —» /q(1710)7 —» KK 7 ) = ( 8 . 5 ^ 9 ) x 10 -4 and 
B(i/)(J/i/)) —> /o(1710)7 —>■ ujiDj) = (3.1 ± 1.0) x 10 -4 . On the other hand, the fraction 
B(ij>(J/ij)) — * /o(171O)7)//5(V , («//'0) -> /o(1500)7) is around order of 10, and /o(1370) is still 
not observed in the J'ljj radiative decays. The mixing matrix can be expressed as m 


/ /o(1370) \ 


/ 0.819(89) 0.290(91) -0.495(118) \ 


( \ nn ) \ 

/o(1500) 

= 

-0.399(113) 0.908(37) -0.128(52) 


\ss) 

V /o(ino) j 


\ 0.413(87) 0.302(52) 0.859(54) j 


\ 1 G) ) 


One can easily see that / o (1370) has a large possibility of the nn component, while / 0 (1500) 
is dominated by the ss component. 

We assume that Glueball component dominates the contribution in V —> /g + 7 since the 
processes from qq components are suppressed by the strong coupling squared a^. Here we 
simply generalize it to all the three scalar mesons. The related branching ratios are given in 
Tab. [TJ where one can see that our results are comparable with data from PDG2014 except 
the predictions for /o(1370). Employing the mixing matrix elements of Eq. (37) based on 
Lattice QCD, we predict a large branching ratio for V(l ) —>■ /o(1370) + 7 while there is 
no signal at experiment. We conclude that the first line in the matrix of Eq. (37) may be 
not precise enough and it need to be checked by the following experiment. 
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TABLE I: The branching ratios (10 4 ) of V(1 ) —>• /q ± 7 and V(1 ) -+ G(0 ++ ) + 7 , where 

V denotes one of heavy quarkonium J/-0, ip(2S) and T (nS), ,/q denotes one of /o(1370), /o(1500) 
and /o(1710). 


Branching ratio (10 4 ) 

This work LQCD ||Tj He et al. [HD] Cheng et al. [3B] PDG2014 [J2] 

B( J/V> ->• /o(1370) + 7 ) a 

9.3 ±2.2 

- 

- 

- 

B{J/ij} -+ /o(1500) + 7 ) 

0.62 ±0.15 

- 

2.9 

1.01 ±0.32 

B(JAM/o(1710) + 7 ) 

28.0 ± 6.6 

- 

14.5 

> 15.6 

B{ J/i/) -+ G( 0 ++) + 7 ) 

38 ± 9 b 38 ± 9 

- 

- 

- 

BM2S)->/ 0 (1370) + 7 ) 

1 47 + 0 - 83 
i - 4 O -0.34 

- 

- 

- 

B{i>(2S) -± / 0 (1500) + 7 ) 

o.97i8:“ 

- 

- 

- 

0^(25)^/0(1710)+ 7 ) 

4-4±?:B 

- 

- 

> 0.9 

B(ip(2S) -+ G(0 ++ ) ± 7 ) 

5-9±?i 

- 

- 

- 

0(T^/o(137O) + 7 ) 

0 q 2 +0 ' 18 - 

u - oz -0.08 

4.8 

- 

- 

B {T -+ / 0 (1500) + 7 ) 

0 021 +0 012 

U.UZi_o 005 

4.2 

- 

< 0.15 

B(T->/o(1710)+ 7 ) 

0 Qfi + 0 - 55 

U.JD_o.23 

1.5 

- 

< 2.6 

B(T -+ G(0++) + 7 ) 

1 3 +0 - 7 

i -°-0.3 

- 

- 

< 2.6 

B(T(2S) —> / 0 (1370) + 7 ) 

0 26 +0 ' 14 - 
u - zd ^0.06 

- 

- 

- 

B(T(2S) -± / 0 (1500) + 7 ) 

0 nifi+ 0 009 - 

U.UlO_ 0 004 

- 

- 

- 

B(T(2S) -+ / 0 (1710) + 7 ) 

n 77 + 0.44 _ 

u - ‘ ‘- 0 .I 8 

- 

- 

< 5.9 

B(T(2S) -+ G(0++) + 7 ) 

1 n+ 0.6 _ 

i - u -0.2 

- 

- 

< 5.9 

B(T(3S) -± / 0 (1370) + 7 ) 

0 30 +0 ' 17 - 

u " 5U -0.07 

- 

- 

- 

B(T(3S) / 0 (1500) + 7 ) 

0 019 + 0010 - 

u.uiy_ 0 00 5 

- 

- 

- 

B(T(3S) ->■ / 0 (1710) + 7 ) 

0 90 +0 ' 51 
u.yu_ 0 .2i 

- 

- 

- 

B(T(3S) G(0++) + 7 ) 

1 2 +a7 

1-Z —0.3 

- 

- 

- 


^Here we adopt the mixing matrix as Eq. (371 from Lattice QCD, however the mixing matrix elements for 


/o(1370) we think still need to be tested further. 
b We use the Lattice QCD result to extract the decay constant of Glueball. 


V. CONCLUSION 

In this paper, we have established the factorization formulae for heavy vector quarko¬ 
nium radiative decays into scalar Glueball, by studying one-loop corrections to the hard 
kernel, LDMEs of heavy quarkonium, and LCDA of scalar Glueball. The NRQCD(^)LCDA 
factorization formulae shall be valid to all orders of the strong coupling constant a s in the 
leading-order twist and heavy quark velocity, after considering a two-dimensional LCDA of 
Glueball which is defined by non-local light-cone gauge-invariant operators matrix elements. 
The universality of LCDA for Glueball ensures us to extract its decay constant. Match¬ 
ing to the CLQCD results, we have extracted the decay constant for scalar Glueball, i.e. 
f g = 0.0386^o;Q04gGeV. We also predict the branching ratios of scalar Glueball from other 
heavy vector quarkonia such as ip(2S) and T (nS), which can be checked in the upcom¬ 
ing experiment. The factorization formulae can also be applied to pseudoscalar and tensor 
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Glueball production from heavy vector quarkonia. A systematic study on Glueball produc¬ 
tion and decay with different quantum numbers shall be investigated in order to hunting 
and identifying Glueball with a large confidence. We will address these issues in following 
studies. 
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Appendix 


In the appendix, we give the explicit results of the short-distance coefficients for hard 
kernels. 

Tja _ i( 2u - l)-Nclog(u) /X 2 (2 u - 1)N c (U 2 (1 - 2 u) + log( m)( log( m) - 2 + 21og(2))) 

- I ~ _ x log 


6 (m — 1 ) 


mi. 


6 (m — 1) 


< 2 “ - < 54 - ^ - 18 - 2/ “> ^ 1 „ ((2u - 1 )(3fl, - 2B 2 + SC, 


144(m — 1 )u 8 (u — 1 )u 

-( 8 u 2 - 8 u + 3)C 2 ) + 2 ( 8 u 2 -9 u + 3 )B 3 )} - u (1 - u ), 

Tja _ Tja 

^1 ~ ^0 ? 


(38) 


m = 


r (2u 2 — 2u + 1) N c log(u) ti 2 (2m 2 — 2u + 1) N c 

[-- - log ~ 2 ~ - 77 —i——-(Li 2 (l - 2 u) + log(rt) (log(rt) - 2 


u 


mf. 


+ 2,o g (2))-#) + YfYW + 


u — 
b\B 2 


+ 


b b 3 B 3 


+ 


b b 4 B 4 


12 u 96(m — 1) 2 m 2 96(m — l)u 48(w — 1 )m 2 (2m — 1) 12(w — l) 2 

c\C i c b 2 C 2 c b C 3 c b CU 


+ 


+ 




192(1 - 2m) 4 (m - 1) 2 m 2 48(m-1) 2 m 2 32(1 - 2m) 4 (m - 1)m 2 24(m - 1) 2 m 

C 5 G 5 c qCq (m 4 — 2 u 3 + u 2 — 3) Cf 


+ 


96(1 - 2w)% 2 48(1 - 2 uY(u - l) 2 

(2m 3 — m 2 — 9m + 1) Cp 


2 (m - 1 ) 2 m 2 


+ 

-c 7 


-4C/4 


4(m — 1)m 2 


-Cg + M — > (1 — m) , 


(39) 
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where 


fi = 12(Li 2 (l - 2 u) + Li 2 (2 u) + log(l - 2w)(log(w) + log(2)) - 7r 2 , 

f\ = ((438 - 4v f 2 )m 4 + (8vr 2 - 876 )m 3 + (585 - 6t t 2 )u 2 + (2vr 2 - 147 )u - 6)N C 

-12(52m 4 - 104m 3 + 46m 2 + 6m - 3 )C F - 90m 4 + 180m 3 - 433m 2 + 343m - 36, 
b\ = 3(180m 4 - 360m 3 + 237m 2 - 57m - 2)N C + 12(2m 4 - 4m 3 + 8m 2 - 6m + 3 )C F 
-1172m 4 + 2344m 3 - 1587m 2 + 415m - 36 , 
b b 2 = —15(1 — 2u) 2 N c + 12C f — 124m 2 + 124m — 45 , 

b\ = 3(1 - 2m) 2 (48m 3 - 92m 2 + 63m - 10) N c + 12(6m 4 - 37m 3 + 47m 2 - 21m + 3 )C F 
+1312m 4 - 2588m 3 + 1740m 2 - 455m + 36 , 

b\ = (6m 2 — 6m + 9) 1V c — 3 (4m 2 + 3m — l) C F — 2u 2 + 3m — 1, (40) 

c\ = 3 (1792m 8 - 7168m 7 + 12688m 6 - 12976m 5 + 8368m 4 - 3472m 3 + 895m 2 - 127m + 8) N c 

+4(m - 1)m(1 - 2m) 4 (51CV - 40m 2 + 40m - l) , 
c\ = -3 (16m 4 - 32m 3 + 27m 2 - 11m + 2) (1 - 2u) 2 N c + 6(1 - 2u) 2 C F 
+m (100m 3 - 200m 2 + 117m - 17) , 

c b 3 = u (-128m 7 + 184m 6 + 932m 5 - 2372m 4 + 2224m 3 - 1039m 2 + 246m - 23) N c 
—4(1 — 2m) 4 (3m 3 — 4m 2 — 5m — l) Cp , 

4 = 3 (10m 3 - 24m 2 + 8m + 7) C F + u (-3 (m 2 - 3m + 4) N c - 32m 3 + 37m 2 + 4m - 9) , 

4 = 3 (128m 7 - 648m 6 + 1356m 5 - 1256m 4 + 472m 3 - 3m 2 - 40m + 8) N c 

+12 (3m 2 - 20m + 24) (1 - 2m) 4 Cf - 4m(9m + 2)(1 - 2m) 4 , 

4 = (u- 1) (3 (128m 6 - 56m 5 - 220m 4 + 232m 3 - 72m 2 + 7m + l) 1V C + 16 (m - 1)(1 - 2m) 4 ) 
-12(1 - 2m) 4 (3m 2 + 7m + 16) C F , 

(41) 


and 

H\ 


H b + 


c\C\ 


+ 


4c 3 


+ 


c\Ca 


+ 


%C 6 


96(1 — 2m) 4 (m — 1)m 96(1 — 2m) 4 (m — 1)m 24(m — 1)m 48(1 — 2 m) 4 m 2 

CgC7 (2 m 4 — 4m 3 + m 2 + m — 3) C F 


24(1 — 2m) 4 (m — l) 2 
(4m 2 — 3m + 5) C F 
2 (m — 1)m 


(m — 1) 2 m 2 
C$ + m y (1 m)] , 


-a 


(42) 
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c\ = -3 (352m 6 - 1056m 5 + 1752m 4 - 1744m 3 + 992m 1 2 3 4 5 6 7 - 296m + 4l) N c 
+24 (4m 2 - 4m + 3) (1 - 2m) 4 C f + 4 (16m 2 - 16m + 45) (1 - 2m) 4 , 
c b 3 = 12(1 - 2m) 4 (2m 2 + 3m + 1) C F - (m - 1)(3(176m 6 - 496m 5 + 796m 4 
-792m 3 + 456m 2 - 129m + 15)N C - 4(1 - 2m) 4 (8m 2 - 5m - 3)), 
c\ = u (—3(4m - 5)N C + 32m 2 - 28m - 45) - 12 (2m 3 - 5m 2 + 3 )C F , 
c\ = m(3 (176m 6 - 560m 5 + 988m 4 - 1016m 3 + 584m 2 - 183m + 28) N c 

-2(1 - 2m) 4 (16m 2 + 10m - 55)) + 12 (4m 2 - 15m + 12) (1 - 2m) 4 C f , 

Cg = (m — 1)(3 (176m 6 - 496m 5 + 796m 4 - 792m 3 + 456m 2 - 129m + 15) N c 

-4(1 - 2m) 4 (8m 2 - 5m - 3)) - 72(1 - 2m) 4 (m 2 + 1 )C F . (43) 

The coefficients 6*, c t and q are related to the scalar Passarino-Veltman integrals defined in 
Ref. [48) 09], and here we have the relation C, = m 2 C { - : 

B 1 = B 0 (0,m 2 ,m 2 ) , 

B 2 = B 0 (-m 2 , 0, m 2 c ) , 

B 3 = B 0 ((1 - 2M)m 2 ,0, m 2 ) , 

R 4 = B 0 (4Mm 2 ,m 2 ,m 2 ) , 

C? = C 0 (-m 2 ,m 2 ,0,m 2 ,0,m 2 ) , 

C 2 — C 0 (0, (1 — 2m)tm 2 , (2m — l)m 2 , m 2 , m 2 , 0) , 

C 6 = C 0 (m 2 ,4Mm 2 , (2 m — l)m 2 ,0,m 2 ,m 2 ) , 

Cl = C 0 (m 2 ,0, (1 - 2M)m 2 ,0,m 2 ,m 2 ) , 

= C 0 (-m 2 ,0,(l-2M)m 2 ,0,m 2 ,m 2 ), 

Cg = C 0 (0,0,4Mm 2 ,m 2 ,m 2 ,m 2 ) , 

Cl = C 0 (4ml, 0,0, ml, ml, ml) , 

Cg = C 0 (4ml, 0,4uml, ml, ml, ml) . 

(44) 
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